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The magnetization plateaux in two dimensionally coupled 5 = 1/2 dimerized zigzag Heisen- 
berg chains are investigated by means of the bond operator mean field approximation. In the 
absence of the interchain coupling, this model is known to have a plateau at half of the satura- 
tion magnetization accompanied by the spontaneous translational symmetry breakdown. The 
parameter regime in which the plateau appears is reproduced well within the present approxi- 
mation. In the presence of the interchain coupling, this plateau is shown to be suppressed. This 
result is also supported by the numerical diagonalization calculation. 
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1. Introduction 

The phenomenon of the magnetization plateau in 
quantum magnets has been attracting broad interest 
as an essentially macroscopic quantum phenomenon in 
which macroscopic magnetization M is quantized to the 
fractional values of the saturated magnetization Mg}^"^^ 
Among them, the plateau associated with the sponta- 
neous translational symmetry breakdown (STSB) is most 
interesting because such a plateau is the manifestation 
of the many body effect.^ 

Theoretically, the plateaux in one dimensional quan- 
tum magnets are best understood, because the size prob- 
lem inherent to numerical studies are less serious and the 
analytical tools such as bosonization and conformal field 
theory are also available. On the other hand, in real ma- 
terials, the interchain coupling is inevitable. ^"^^'^ How- 
ever, the limitation of the system size becomes serious 
in the numerical studies of the systems in higher dimen- 
sions. Therefore it is desirable to to develop a reliable 
approximation scheme appropriate for the description of 
the plateau state in higher dimensions. 

Among various approximation schemes which are com- 
monly used in quantum spin systems, the bond oper- 
ator mean field approximation (BOMFA) is known to 
be powerful for the description of the spin gap state in 
the absence of magnetization. ^^^^^ Considering that the 
plateau state can be regarded as a spin gap state with 
non-zero magnetization, this method must be also useful 
for the description of the plateau state. Actually, there 
are several attempts to describe the magnetized states in 
terms of the bond operators.^' ^^"^^ In the present work, 
we propose the application of the BOMFA to the plateau 
state in the two dimensionally coupled S — \/2 dimer- 
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ized zigzag Heisenberg chains. Especially, we show that 
the plateau accompanied by the STSB can be described 
by the BOMFA with corresponding configuration of sin- 
glet and triplet pairs. 

The single S = 1/2 dimerized zigzag Heisenberg 
chain is known to have a magnetization plateau at m(= 
M/Ms) = 1/2 in the appropriate parameter range. 
Our approximation quantitatively reproduce the phase 
diagram at m = 1/2. Motivated by the success in one 
dimension, we further investigate the effect of the in- 
terchain coupling within our approximation scheme. It 
is found that the interchain coupling suppresses this 
plateau. The numerical diagonalization results also sup- 
ports this conclusion. 

This paper is organized as follows. In the next sec- 
tion, we briefly review the bond operator transforma- 
tion. The Hamiltonian of the two dimensionally cou- 
pled S — 1/2 dimerized zigzag Heisenberg chains is pre- 
sented and rewritten using the bond operators in §3. 
In section 4, the mean field approximation appropriate 
for the plateau state is introduced and applied for the 
present Hamiltonian. The condition for the stability of 
the plateau is derived. The ground state phase diagram 
is presented in §5. The last section is devoted to summary 
and discussion. 

2. Bond Operators 

The bond operators s, s^, ta and — ±,0) are 

defined for a pair of two 5 = 1/2 spins 5'l and 5'r as^^"^^ 

^(in>-UT», 

= - I TT) , (2) 

= IU>, (3) 
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where |0) is the vacuum state and \SIS^) is the eigen- 
state of SI and 5^. These operators satisfy the bosonic 
commutation relations, 



[s,st] =1, 



Oct J 1^ r. 



Saf}, [s,ti] = 0. 



(5) 



Because the physical state of each spin pair is either sin- 
glet or one of the three triplet states, these boson oper- 
ators must satisfy the following condition, 

(6) 



sh + tlt+ + tlt_ + 4*0 = 1- 



The components of the spin operators S'r and S'l are 
rewritten in terms of the bond operators as follows: 

= ^ {± i^^t- + t+^s) - {tlto - 4i-) } , (7) 

Sr,l = 71 {± (^^*+ + *-^^) + - 4*+) } , (8) 

^R,L = ^ {± (s^to + toh) + (tlt+ - tit.) } , (9) 

where upper and lower sign in rhs correspond to R and 
L in Ihs, respectively. 

3. Hamiltonian 

The Hamiltonian of the two dimensionally coupled S = 
1/2 dimerized zigzag chains in the magnetic field H along 
^;-direction is given by, 

N 

W = J^^{(l + (5) S^{t,l)-Sn{t,l) 

1=1 i=l 

+ {l-S)Sj^{i,l)-Sn{i + l,l) 

+ Ai (SL(i, I) ■ Si^ii + 1,0 + Sn{i, I) ■ Snii + 1, /)) 

+ X2{SR{i,l)-SL{i,l + l) 



+ Snii,l)-Sj^{i + 1,1 + 1))} 

Na N 

1=1 i=l 



(10) 



where N is the length of each chain and Nc is the num- 
ber of chains. The chains are distinguished by the index 
I. The dimers coupled via the stronger nearest neigh- 
bour bonds J(l + S) are the elementary components of 
the present model. In the following, these spin pairs are 
simply called 'dimers'. They are distinguished by the in- 
dex i in each chain. The indices R and L distinguish the 
two spins in each dimer. Accordingly, the spin operators 
SB.{i,l) and Si,{i,l) are the right and left spins in the 
i-th dimer on the l-th chain, respectively. The weaker in- 
trachain nearest neighbour coupling, the intrachain next 
nearest neighbour coupling and the interchain coupling 
are denoted by J(l — S), JAi and JA2, respectively, as 
depicted in Fig. 1. This model can be also regarded as 
a dimerized triangular Heisenberg model. The uniform 
triangular lattice corresponds to Ai = A2 = 1 and 6 = 0. 

The Hamiltonian (10) can be represented by the bond 
operators for the dimers connected by J(l + 6) bonds as 




J(1+5) 
J(1-5) 



Fig. 1. Two dimensionally coupled dimerized zigzag chains. The 
filled circles represent the 5 = 1/2 spins. 



follows. 
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+ ^ E ti{i,l)ta{i,l)]-Ki,l)(s^{i,l)s{i,l) 
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+ E 4(^,0*a(^,0-l)l, (12) 

= - E E I (1 - 0*o(i, + l)s{i, I) 

1=1 i=l 

+ tl{i, l)t+{i, I) - tl{i, l)t-{i, I)) 

X (^s\i + 1, l)to{i + 1,0+ 4(* + ^,l)s{i + 1, 

- tl{i + 1, l)t+{i + 1,0+ i-(« + 1, l)t-{i + 1, 0) 
+ \^(^s''ii,l)t-{i,l) + tl{i,l)s{i,l) 
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= ^ ^ -Ai l)to{i, I) + tlii, l)s{i, 0) 

1=1 i=l 

X (^s\i + 1, l)to{i + 1,0+ 4(« + 1' Os(« + 1, 0) 

+ (tl{i,l)t+{i,l)-tl{i,l)t-{i,l)') 

X (tl{i + 1, Ot+(i + 1, /) - t-{i + 1> 0*-(« + 1. 0) 

+ {(st(z,;)i_(z,0 + tU^'Os(*,0) 

X (st(j + 1^ + 1, Z) + tl{i + 1, Os(i + 1, 0) 
+ (tl{i,l)t-{i,l) - tl{i,l)toii,l)^ 

X (4(» + 1' Oi+(« + 1,0- + 1> Oio(« + 1, 0) 
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= - ^ ^ ^A2 Oto(i, + 4(*> Os(*, 

i=l i=l 

+ tl{i,l)t+{i,l)-tl{i,l)t-{i,l) 

+ s^(i + 1, /)to(i +1,0 + + 1, + 1, 

+ tl{i + l,l)t+{i + l,l)- tl{i + l,l)t_{i + l,l)^ 

X (st(j^ I + / + 1) + 4(i, / + l)s{i, I + 1) 

- tlii, I + l)t+(i, l + + tl{i, I + l)t-{i, I + 1)) 
+ ^^(^s''{i,l)t-{i,l) + tl{i,l)s{i,l) 

- tl{i,l)t-{i,l) + tl{i,l)to{i,l) 

+ s^i + 1, l)t-{i + 1,0 + + 1, Os(« + 1, 

- tlii + 1, Ot-(« + 1,0 + + 1> 0^o(« + 1, 0) 

X (fit / + l)t+{i, l + tl {i, I + l)s{i, I + 1) 

- tl{i, I + l)t+{i, l + + tl{i, I + l)toii, I + 1)) 
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(b) 



Hi = (t;(i,Ot+(z,0 - iL(z,Oi-(^,0) , (16) 

where the Lagrange multipliers /z(i, are introduced to 

account for the constraint (6). 

4. Bond Operator Mean Field Approximation in 
the Plateau State 

On the m = 1/2 plateau state, it is expected that half 

of the dimcrs on the J(l + (5)-bonds are in the singlet 
states and others are in the triplet states polarized along 
the ^-direction as far as 1 + 5 >> 1 — 5, Ai, A2. In this 
state, the translational symmetry i — > i + 1 is sponta- 
neously broken to i ^ i + 2. In the two dimensional case. 



Fig. 2. Two possible configurations of the triplet and singlet pairs 
in the m = 1/2 plateau state. The spin pairs in each oval form 
singlet pairs and the spins represented by the thick up arrows 
are spins in the triplet state polarized along the magnetic field. 



there arc two possible singlet-triplet configurations as de- 
picted in Fig. 2(a) and (b). We have carried out our cal- 
culation for both configurations. However, it turned out 
that the configuration (a) gives wider plateau. Therefore, 
we only present the calculation for the configuration (a) 
in the following. The calculation for the configuration (b) 
can be carried out almost in the same manner. 

Based on the above explained physical picture of the 
plateau state, we assume the condensation of s{2i, I) and 
t-i-(2i+l, bosons and keep the ground state expectation 
values of these operators as s =< s^(2i, >=< s(2i, > 
and t =< t\{2i + l,l) >=< t+{2i + l,l) >. 

Wc also make global approximation for the Lagrange 
multiplier iJ,{i,l). In the present case, the (2i,0-th sites 
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and the (2z + l,Z)-th sites are inequivalent. Therefore + 
we assume different values /i° and /i° for fj,{2i, I) and 
/x(2i + 1,1), respectively. Thus, the ground state energy ^j^gj-e 
is obtained as, 



Eq 



NrN 



(17) 



neglecting all quantum fluctuation terms. 

The parameters s, t, ji'^ and /U° are determined so as 
to optimize Eq as follows. 



J{l + 5), 



(18) 



4 

t = 1 



Neglecting the terms higher than the third order in 
t±{2i, I), to{2i, I), s(2i + 1, 1), to{2i + 1, and t-{2i + l, I) 
and carrying out the Fourier transformation. 



ta{2i,l) 



N,N 



ik-r{2i,l) 



, (a = 0,±), 
(21) 



s{2i + l,m) = 
i„(2i + l,0 = 



NrN 



Es°(fe)e^'=-^(2i+i,;)^ (22) 



NoN 



ife-T-(2j+l,i) 



(a = 0,-), (23) 



where r{i, I) is the position of the {i, Z)-th dimer, we ob- 
tain the BOMFA Hamiltonian as, 



n 



MF 



Eq 



+ A3(fe)t^ot(fe)ie(fe) + ^°(fe)s;t(fe)gO (fe) 

+ A° (fe)<°-f(fc)i° (fe) + A°(fc)i°^(fc)i°(fc)) 

+ ^Pik){s°Hk)(f^{-k)+tt{k)) 

k 

s°{k) (ti{-k) + f_!ik))} 

^iQ{k){t°\k) (t^_^{-k) + tUk)) 

k 

m (t;(-fc)+d(fe))} 
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l{ti{kn{-k)+ti\k)ti\-k))^ 
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A^(fc) = J(l + 5)-i/ + JAi + ^(l-(5) 



+ ^(l-cosfcy). 



^l(fe) 



J, 



(19) Ag(fc) 
(20) 



j(i + 5) + ir- jAi- -(i-(5) 

^(1 + COSfcy), 

J(l + 5) - COSfcy, 

JAo 



A°(fe) = -J{l + 6) + H - ^{1 + cosky), 
A°{k) = H cos ky), 



A'Lik) = 

P{k) = 

Q{k) = 

S{k) = 



2H-JX2, 

'1-S 



Ai I Jcosfcj; - cos(fca, - ky), 



^—^Jsinkx - ■^^sin(fcx - ky), 
^^^sinfcy, R{k) = cos ky. 



where k^ and ky are the momentum components conju- 
gate to i and I. 

In the matrix representation, eq. (24) is rewritten as. 



H 



MF 



EQ + Y^u\k)D{k)u{k) 

k 

k 

+ Al{k) + A'^_{k) + Al{k)), 

where 

M^(fc) = ( t\k) t{-k) ) 

with 
<t(fc) = 

( fl{k) tl\k) fl{k) s°^k) t°\k) t°l{k) 



(25) 
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D'i(fc) = 

/ A%{k) 

Alik) 
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(28) 



-(^) - ( ) • (-) 

The summation is taken over a half of the Brillouin 

k 

zone (say > 0).lt is known that this type of Hamilto- 
nian is diagonahzable by the Bogohubov transformation 
as, 

Y^t\k)D{k)t{k) ^M(fe) (alik)a^{k) + 

fe k (U=l,6 ^ ^ 

(30) 

with six modes a^(fe), {n = 1, ..,6) which have positive 
excitation energies huj^{k) if and only if the matrix D{k) 
is positive definite for all allowed values of fc.^^ From this 
condition, the plateau region is determined numerically. 

5. Results 

5.1 Single Dimerized Zigzag Chain 

Figure 3 shows the magnetic phase diagram of the 
S = 1/2 one-dimensional dimerized zigzag chain (A2 = 0) 
on the H — 5 plane determined by the present method 
for Ai = 0.2. Compared to the numerical results (Fig. 
3 of ref. 4), the BOMFA overestimates the width of the 
plateau. Especially, the plateau width remains finite even 
for S = where no plateau is observed in numerical cal- 
culation. Nevertheless, the approximate position of the 
plateau and the range of 6 for which the plateau appears 
is fairly well reproduced. 

The latter feature is visible more clearly in Fig. 4 which 



shows the ground state phase diagram on the 5 — Xi 
plane. The isolated dimer limit, which is the starting 
point of the BOMFA, is characterized by Ai = and 
= 1.0. For finite Ai and/or 1 — ^, the interdimer in- 
teraction is switched on and the many body effect stabi- 
lizes the TO = 1/2 plateau state accompanied by STSB. 
Comparing the present results with the numerical phase 
diagram by Tonegawa and coworkers,'^' ^(Fig. 1 of ref. 4) 
the plateau region is well reproduced by the BOMFA. 
This implies that the present approximation takes into 
account the many body effect appropriately. 

5.2 Two Dimensional Coupled Dimerized Zigzag Chain 
The results of the preceding subsection indicates that 
the BOMFA gives reliable results even for the one- 
dimensional case in which the mean field type approx- 
imation is expected to be rather poor. Encouraged by 
this success, we further investigate the effect of the in- 
terchain coupling in this subsection. 

The TO = 1/2 ground state phase diagram on the 
S — Xi plane is shown for various values of A2 in Fig. 
5. It is seen that the interchain interaction suppresses 
the TO = 1/2 plateau. Actually, no plateau appears for 
A2 > 0.634. This conclusion is physically reasonable be- 
cause the present plateau is understood as the solidifi- 
cation of the singlet and triplet pairs which alternate 
on the strong bonds. ^ The interchain coupling tends to 
align the spins and destroy the singlet pair formation. 
Therefore the present plateau is destabilized by the in- 
terchain interaction. Considering that the mean field ap- 
proximation becomes more reliable as the dimensionality 
becomes higher, we expect the present calculation is re- 
liable in the coupled chain system. 

To confirm the above conclusion numerically, we have 
carried out the numerical diagonalization calculation for 
Nc = N ~ 4 and 6 clusters. Figure 6 shows the Nc)- 
dependence of the magnetic fields at the lower end Hd 
and upper end Hc2 of the m = 1/2 plateau for (a) 
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Fig. 5. The ground state phase diagram of the two dimensionally 
coupled 5 = 1/2 dimerized zigzag Heisenberg chains for m = 
1/2. The region encircled by each curve and Ai-axis is the plateau 
region. The values of A2 are indicated for each line. 



,2 




(a) 5=0.5 ^1=0.4 
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(b) 5=0.6 
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^2=0.1 
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A,2=0.4 


♦ 



0.1 



0.2 
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Fig. 6. The magnetic fields at lower end Hd and upper end Hrjz 
of the m = 1/2 plateau of the 5 = 1/2 coupled dimerized zigzag 
Heisenberg chains for the clusters with A'^ = A^c = 4 and 6 plotted 
against l/N for (a) (<5,Ai) = (0.5,0.4) and (b) (0.6, 0.5). The 
solid line is the linear extrapolation in 1/N. The magnetic field 
is scaled by J. 



{5,\i) = (0.5,0.4) and (b) ((5,Ai) = (0.6,0.5), respec- 
tively. Of course, the system size dependence of Hd and 
i?c2 is not clear, so that there is large ambiguity in the 
extrapolation to iV ^ 00 from only two values of N . 
Therefore the absolute values of Hd and Hc2 are not 
reliable. However, it is clear from Fig. 6 that the extrap- 
olated values of the plateau width decrease with A2 . It is 
also checked that this tendency is insensitive to the way 
of extrapolation. Therefore the numerical diagonalization 
results support the validity of the BOMFA result. 

6. Summary and Discussion 

In the present work, we have investigated the m = 1/2 
plateau of the coupled dimerized zigzag chain by apply- 
ing the BOMFA to the plateau state. It is found that 
the m = 1/2 plateau associated with STSB can be de- 
scribed by BOMFA. In the single dimerized zigzag chain, 
the ground state phase diagram obtained by the numer- 
ical calculation is reproduced quantitatively well. In the 
coupled dimerized zigzag chains, the plateau is found to 
be suppressed by the interchain coupling. This is also 
consistent with the numerical diagonalization results. 

Similar conclusion is also obtained by Kolezhuk^ sev- 
eral years ago. However, his calculation treats the in- 
terdimer interaction only perturbatively and no quanti- 
tative phase diagram is obtained. Furthermore, the ef- 
fective kinetic energy of triplet pairs are assumed to be 
much smaller than the interaction energy between these 
pairs, although both energies are in fact the same order 
in Ai, A2 and 1 — 5. In the present treatment, the both 
energies are equally taken into account. 

Sommer and coworkers^^ and Matsumoto and cowork- 
ers^' also applied BOMFA to the magnetization process 
of the quantum magnets. However, these authors con- 
centrated on the non-plateau part of the magnetization 
curve and explained the field induced transverse antifer- 
romagnctic ordering in two and three dimensional mod- 
els. The simple-minded application of their method to 
the one-dimensional case also leads to the transverse or- 
dering. This is obviously unphysical, because the ground 
state of the quantum spin chains in the non-plateau re- 
gion should be the Luttingcr liquid state. Furthermore, 
the plateau associated with the STSB cannot be de- 
scribed by their scheme, because STSB occurs only on 
the magnetization plateau. On the other hand, our ap- 
proach concentrate on the plateau state and is not suit- 
able for the description of the non-plateau part of the 
magnetization curve. It is, however, applicable to the 
one-dimensional case and can describe the plateau state 
with STSB. Therefore we may conclude that our method 
and the method of refs. 9, 18, 19 are complementary to 
each other. 

For (5 = and A2 = 0, Okunishi and Tonegawa^'^ 
have shown that the m = 1/3 plateau appears in the 
strongly frustrated regime. For (5 = 0,Ai = A2 = 1, 
the present model reduces to the triangular lattice for 
which the rn = 1/3 plateau is believed to be present. ^^'^^ 
In this context, it is interesting to investigate the inter- 
play of two dimensionality, dimerization and frustration 
in the m = 1/3 plateau problem. We have also tried to 
apply the present method to the m = 1/3 plateau. How- 
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ever, the width of the plateau region on the 5 — Ai plane 
determined by the BOMFA is substantially wider than 
the numerical results by Okunishi and Tonegawa^'* for 
5 = and preliminary DMRG calculation by one of the 
authors (KH) for (5 0. The origin of this discrepancy 
is the following. As discussed by Okunishi and Tone- 
gawa/'^ this m = 1/3 plateau is basically of classical 
origin with up-up-down local spin configuration realized 
in the Ising limit rather than the quantum mechanical 
plateau resulting from the singlet pair formation on the 
strong bonds. Therefore the present approach based on 
the singlet pair formation is not suitable for the descrip- 
tion of the m = 1/3 plateau. The same is true for the 
triangular lattice {6 = 0, Ai = A2 = 1) which also has 
the m = 1/3 plateau. Therefore in order to apply 
the present method to the m = 1/3 plateau, it would 
be necessary to refine the approximation so as to incor- 
porate the effect of the ground state polarization of the 
singlet sites. This is left for future studies. 

In conclusion, the BOMFA is a promising tool for the 
investigation of the magnetization plateau of quantum 
origin including those induced by many body effect asso- 
ciated with the translational symmetry breakdown. The 
application of the present method to wide range of mod- 
els is hoped in the future. 
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